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ABSTRACT 

^^The  acoustic  pressure  in  vater  due  to  a plane  wave 
incident  on  an  elastic  plate  is  considered.  The  plate  is  backed 
by  air  and  by  an  infinite  array  of  parallel  stiffeners  (or  mechanical 
resonators).  The  acoustic  pressure  on  the  plate  is  calculated  as  a 
function  of  the  spacing  between  the  stiffeners  or  resonators,  the 
plate  thickness,  the  angle  of  incidence,  and  frequency. 
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INTRODUCTION 

In  1961  TRG  showed'*'  that  the  observed  inability  to  steer 

a planar  sonar  array  to  endfire  is  caused  by  the  loss  in  pressure 

at  the  hull  of  the  ship  due  to  the  finite  rigidity  of  steel  plates 

of  practical  thickness.  This  theoretical  deduction  was  later 

2 1 

supported  by  experimental  measurements.  TRG  suggested  that  the 
use  of  mechanical  resonators  (also  called  stiffeners)  fastened  to 
the  surface  of  the  hull  could  increase  the  impedance  of  the 

hull  to  a satisfactory  level  over  a bandwidth  sufficient  for 
practical  purposes.  The  effectiveness  of  resonators  was  later 
verified  experimentally. ^ 

Theoretical  studies  of  the  effectiveness  of  resonators 

have  considered  either  single  resonators^  or  resonators  distributed 

uniformly  and  continuously^  behind  an  infinite  elastic  plate.  The 

study  reported  here  will  be  concerned  with  an  infinite  set  of 

resonators  or  parallel  stiffeners  distributed  at  regular  intervals 

behind  an  infinite  plate  in  contact  with  water.  The  results  for 

this  discrete  distribution  approach  the  results  for  the  uniform 

and  continuous  distribution  when  the  spacing  between  resonators 

approaches  zero.  The  purpose  of  this  theoretical  study  is  to  find 

the  maximum  spacing  between  resonators  at  which  the  plate  is  still 

sufficiently  rigidized.  We  will  calculate  the  pressure  at  the 

plate  due  to  an  incident  plane  wave;  the  plate  will  be  considered 

to  be  sufficiently  rigidized  when  this  pressure  at  the  plate  does 

not  differ  appreciably  from  the  pressure  at  an  ideally  rigid  plate. 

It  has  been  shown  experimentally  that  a plate  can  be  rigidized 

2 5 

with  practical  spacings  between  resonators;  ’ however,  we  may 
have  used  smaller  spacings  and  more  resonators  than  necessary, 
which  motivates  this  study.  One  purpose  of  this  work  is  to  provide  a 
more  thorough  analytic  foundation  for  the  design  of  resonant 
stiffeners . 

The  numerical  results  are  presented  in  the  next  section. 

The  details  of  the  mathematical  derivations  are  given  in  two 
appendices . 
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NUMERICAL  RESULTS 


Consider  an  elastic  plate  of  thickness  h in  contact 
with  water  occupying  the  half-space  x 0 , see  Figure  1.  Resonators 
or  parallel  stiffeners  are  attached  to  the  plate  at  ns-b  < y < ns , 
i = 0 , +1 , +2  where  b is  the  width  of  the  resonator  attachment 

(assumed  to  be  very  small  as  compared  to  a wavelength  in  water  or 
flexural  wavelength  in  the  plate) , and  s is  the  spacing  between 
resonators.  The  resonators  are  assumed  to  be  infinite  in  the 
z-direction  (normal  to  the  plane  of  the  paper  in  Fig.  1). 

The  equations  describing  the  pressure  at  the  surface 
of  the  plate  and  the  plate  velocity  in  the  x-direction  are  derived 
in  Appendices  A and  B (both  appendices  deal  with  the  same  problem, 
but  from  slightly  different  points  of  view).  The  effects  of  reson- 
ators are  given  by  specifying:  1)  the  force  which  the  resonator  ex- 
erts on  the  plate  at  the  attachment  point;  this  force  is  proportion- 
al to  the  plate  displacement;  2)  the  moment  which  the  resonator  ex- 
erts on  the  plate;  this  moment  is  proportional  to  the  slope  of  the 
plate  at  the  resonator  attachment  point.  Various  designs  of  reson- 
ators are  possible;  for  example,  a simple  and  effective  device  used 
frequently  consists  of  a short  stud,  which  acts  as  a spring,  with  a 
larger  piece,  which  acts  as  a mass,  at  the  end.  An  infinite  number 
of  different  resonator  designs  will  give  the  same  force  and  moment 
at  the  attachment  point;  therefore  for  simplicity  in  the  calcula- 
tions presented  below  we  assumed  that  the  resonator  forces  and  mo- 
ments are  1/10  of  those  for  a simple  steel  rod  of  circular  cross 
section  with  the  first  axial  and  second  transverse  vibration  reson- 
ance at  2500  cps.  While  the  dimensions,  which  are  determined  by  the 
resonant  frequency,  of  such  a simple  resonator  then  must  be  very 
large  (length  42  cms,  radius  6 cms) , one  should  remember  that  an 
infinite  number  of  other  designs  will  give  the  same  forces  and 
moments,  and  a more  efficient  resonator  (such  as  the  stud  plus 
mass  mentioned  above)  will  be  much  smaller  in  size^.  Note  also 
that  we  are  actually  considering  a resonator  only  1/10  of  the  size 
of  the  rod,  since  the  forces  and  moments  are  assumed  to  be  1 10  of 


er 


TM-67-8 


those  for  the  rod.  Most  calculations  were  performed  with  the 
1/10  simple  rod  resonator  only  because  the  mathematical  expres- 
sions for  forces  and  moments  were  easy  to  obtain  and  easy  to 
compute . 

Figure  2 shows  the  average  pressure  magnitude  on 
the  plate  vs.  the  spacing  between  resonators.  The  plate  is  1/4 
inch  thick,  and  the  frequency  is  2500  cps . An  approximation 
in  which  the  resonators  are  assumed  to  be  distributed  uniformly 
and  continuously^  behind  the  infinite  plate  gives  0 db  for  the 
pressure  for  all  spacings  and  angles  in  Fig.  2.  The  pressure 
reduction  for  s ;>  7 cms  is  easily  understood:  the  wavelength 
for  free  flexural  vibrations  is  14  cms  (for  f = 2500  cps,  h = 

1/4  inch),  therefore  stiffeners  placed  at  half-wavelength  in- 
tervals (s  = 7 cms)  cannot  be  very  effective.  While  we  have 
plotted  the  average  pressure  in  Fig.  2 (and  also  in  the  sub- 
sequent figures)  instead  of  the  actual  pressure  distribution 
on  the  plate,  the  actual  distribution  rarely  differs  by  more 
than  0.01  db  from  the  average  pressure,  unless  the  pressure 
levels  are  very  low. 

Figure  3 shows  the  average  pressure  vs.  the  angle 
of  incidence  0 for  resonator  spacings  s=4,5,  and  6 cms.  Figures 
2 and  3 indicate  that  a resonator  spacing  of  4 cms  is  accept- 
able; the  pressure  is  reduced  by  less  than  2 db  for  0^2° . 

Figure  4 shows  the  average  pressure  vs.  the  fre- 
quency. The  corresponding  resonator  or  stiffener  strength 
vs.  the  frequency  is  shown  in  Table  I,  where  F and  G are 
quantities  proportional  to  the  force  and  moment  exerted  by 
the  resonator  as  defined  in  Appendix  B,  Eqs.  (B-12),  (B-13), 
(B-35)  and  (B-36).  Although  the  stiffener  itself  exerts  the 
maximum  force  and  moment  on  the  plate  at  2500  cps,  as  shown 
in  Table  I,  at  0 = 5°  Figure  4 shows  that  the  plate  is  most 
rigid  between  2200  cps  and  2400  cps,  and  the  pressure  is  lower 
at  3000  cps  than  at  2000  cps.  Consequently,  an  improvement  in 
the  pressure  level  in  a total  bandwidth  of  1000  cps  could  be 
obtained  by  either  shifting  this  band  to  lower  frequencies 
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(from  1800  cps  to  2800  cps)  with  a stiffener  resonance  at  2500 
cps  as  in  Figure  4,  or  by  changing  the  stiffener  resonance  frequency 
to  2700  cps  and  keeping  the  band  from  2000  cps  to  3000  cps.  In 
either  case  one  would  not  lose  more  than  2 db  in  pressure  at  the 
extreme  frequencies  of  the  band  with  the  stiffener  spacing  s = 4 
cms  and  0 ^ 5°. 

Figure  5 shows  the  average  pressure  vs.  the  thick- 
ness of  the  plate  for  a resonator  spacing  of  4 cms. 

Figure  6 shows  the  average  pressure  vs.  the  resonator 

strength,  where  we  have  fixed  the  relative  values  of  F and  G in 

_ -j 

such  a way  that  G = -2  x 10  F.  For  small  values  of  F and  G the 
pressure  magnitudes  approach  the  values  for  a 1/4  inch  thick  plate 
without  resonators,  which  are: 

-27.0  db  at  0 = 5°, 

-12.1  db  at  0 = 30°, 

-6.8  db  at  0 = 90°, 

For  a 1/2  inch  thick  plate  without  resonators  the 
pressure  levels  are: 

-21.1  db  at  0 = 5°, 

-6.9  db  at  0 = 30°, 

-2.9  db  at  0 = 90°, 

Figures  7-10  show  the  same  data  for  a half  inch  thick  plate  as 
in  the  previous  figures  for  a quarter  inch  thick  plate. 

Since  the  wavelength  for  free  flexural  vibrations  for 
the  1/2  inch  thick  plate  is  20  cms,  the  pressure  vs.  resonator 
spacing  in  Fig.  7 shows  a loss  at  the  half-wavelength  spacing, 
s = 10  cms.  Very  small  losses  are  shown  in  Fig.  8 at  grazing 
angles,  which  should  be  compared  with  Fig.  3 for  the  1/4  inch 
thick  plate.  Figures  7 and  8 indicate  that  at  f = 2500  cps  a 
resonator  spacing  of  6 cms  for  the  1/2  inch  thick  plate  is 
acceptable. 

However,  if  we  now  compare  Figures  9 and  4,  which 
show  pressure  vs.  frequency  at  a resonator  spacing  of  4 cms,  we 
do  not  see  any  improvement  with  the  thicker  plate  If  we  compare 
Figures  10  and  6,  which  show  pressure  vs.  the  resonator  strength, 
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we  see  that  the  thicker  plate  has  a higher  pressure  level  at  low 
resonator  strengths  and  at  high  resonator  strengths,  but  that 
at  intermediate  resonator  strengths,  where  the  transition  from 
a flexible  to  a rigid  plate  takes  place,  the  slope  of  the  pressure 
curve  is  steeper  for  the  thinner  plate,  and  as  one  increases  the 
resonator  strength,  the  thinner  plate  becomes  rigidized  sooner 
than  the  thicker  plate.  Apparently  the  rigidization  of  the  plate 
depends  on  the  ratio  of  the  plate  mass  to  the  resonator  strength. 
This  explains  the  lack  of  improvement  in  Figure  9 as  compared  to 
Figure  4,  because  at  the  extreme  frequencies  of  the  band  the 
resonator  strengths  are  lower  than  in  the  middle  of  the  band,  see 
Table  I,  and  the  thinner  plate  is  actually  better  rigidized  than 
the  thicker  plate. 
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TABLE  I 


of  Resonator  or  Stiffener  Strength  vs.  Frequency 


Frequency 

cps 

F, 

10  2 
10  newtons/m 

G, 

10 ^ newtons 

2000 

1.14 

-0.85 

2100 

1.51 

-1.49 

2200 

2.10 

-2.52 

2300 

3.23 

-4.48 

2400 

6.30 

-9.65 

2500 

46.83 

-64.71 

2600 

-9.49 

19.01 

2700 

-4.48 

9.34 

3800 

-2.99 

6.58 

2900 

-2.27 

5.26 

3000 

-1.84 

4.48 
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CONCLUS IONS 


The  calculations  show  that  resonators  or  stiffeners 
spaced  at  distances  of  4 cms  (or  less)  rigidize  a 1/4  or  1/2 
inch  thick  plate  sufficiently  so  that  the  pressure  levels  are 
acceptable  over  a bandwidth  of  1000  cps. 

One  should  note  that  in  an  actual  baffle  design  the 
spacings  between  resonators  will  probably  be  irregular,  in  which 
case  the  pressure  levels  may  be  improved,  because  a regular 
spacing  (which  was  assumed  in  the  above  calculations)  permits 
the  development  of  flexural  waves  in  the  plate  and  supports 
various  resonance  effects  which  may  be  suppressed  by  an  irregular 
spacing. 

One  might  think  that  the  performance  of  the  plate 
could  be  just  as  well  improved  by  adding  the  mass  of  the  resonators 
to  the  plate,  thus  making  the  plate  thicker.  However,  an  increase 
in  the  thickness  of  the  plate  also  increases  the  flexural  wave- 
length, and  if  the  wavelength  in  water  divided  by  cos  0 matches 
the  flexural  wavelength,  the  pressure  at  the  plate  is  zero. 

For  example,  if  we  had  added  the  resonator  mass  at  4 cm  spacing 
to  the  1/2  inch  plate,  then  at  2500  cps  the  wavelength  in  water 
would  match  the  flexural  wavelength  in  the  plate  at  endfire,  and 
the  thick  plate  would  act  as  a pressure  release  surface  near 
endfire . 
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APPENDIX  A 

Vibrations  of  a Plate  with  Parallel  Stiffeners 

D.  Yarmush 


The  following  discussion  is  largely  abstracted  from 

Ref.  6. 

We  first  derive  the  equation  of  motion  of  a uniform 
infinite  plate  that  bears  a continuous  but  non-uniform  distribution 
of  axial  and  transverse  resonators.  These  are  idealized  as  being 
attached  at  points,  rather  ti.an  over  areas.  The  plate  is  in  the 
(y,z) -plane. 

Let  A(y,z)  be  the  density  of  axial  point  resonators  at 
the  point  (y,z).  Each  resonator  is  taken  as  of  unit  strength. 

This  means  that  the  force  exerted  on  an  element  of  area  dy  dz  in 
response  to  the  displacement  q is  t,  A(y,z)dy  dz.  (It  is  understood 
that  the  time  dependence  has  been  factored  out  of  q,  and  that  A 
depends  on  the  frequency.)  Let  P(y,z)  be  the  net  pressure  of  the 
water  on  the  plate.  Then  the  equation  of  motion  is 

D(V4q  - k4q)  = P(y , z)  + A(y,z)q 

where  D is  the  plate  stiffness,  and  k is  the  wavenumber  for  free 
flexural  vibrations. 

For  transverse  resonators,  the  motions  in  perpendicular 
planes  will  be  handled  independently.  That  is,  each  physical 
resonator,  which  can  vibrate  in  many  planes,  is  considered  as 
resulting  from  the  superposition  of  two  kinds,  say  the  xy  and  xz 
kinds.  An  xy  resonator  responds  only  to  the  slope  of  the  cross- 
section  of  the  plate  that  lies  in  the  xy  plane,  and  the  moment  it 
exerts  on  the  plate  lies  in  the  same  plane.  Similarly  for  the  xz 
resonators.  These  two  types  are  assumed  to  have  the  same  density 
T(y,z)  over  the  plate. 
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We  will  study  the  effect  of  the  xy  resonators  only. 

Then  the  dimensionality  of  the  system  is  lowered,  and  thin-rod 
theory  can  be  used. 

Consider  a beam  with  a continuous  distribution  of  trans- 
verse resonators.  The  distribution  is  described  by  a density 
function  T(y)  whose  units  are  moment/unit  rotation/unit  length. 
On  an  increment  dy  of  the  beam  the  increment  of  moment  dm  is 
therefore : 


dm  = T (y)  dy 


But  from  simple  beam  theory, 

El  = m, 
dy2 

and  therefore  the  contribution  of  the  transverse  resonators  to 
the  shear  force  per  unit  length  is 

ET  H = — A = (t(v)  drx 

LI  dy4  dy  My'  dy  my'  dy;  * 


For  the  xz  resonators,  the  contribution  is,  by  a similar  argument, 


HZ  (T(z)  Hz> 

The  sum  of  the  two  terms  can  be  written  in  the  form 
div  (T(y,z)  grad  rj(y,2>) 

Therefore  the  equation  of  motion  for  a plate  bearing  axial  point 
resonators  with  strength-density  A(y,z)  and  transverse  point 
resonators  with  strength-density  T(y,z)  is 

D(V^-k^)  rj  = P + At,  + div(T  grad  t,)  , 


where  P is  the  driving  pressure  distribution. 


iy 


TM-67-8 


We  now  consider  the  special  case  in  which  the  distribu- 
tions A and  T are  concentrated  along  the  lines  y = ns , for 
n = 0,  +1,  +2,  and  have  uniform  linear  density  along  these 

lines.  The  distance  s can  then  be  called  the  resonator  spacing. 

The  axial  strength  per  unit  length  at  frequency  to  will 
be  F(o>)  and  the  transverse  strength  G(co)  . We  also  assume  that 
the  pressure  of  the  water  on  the  plate  is  constant  along  each 
line  y = const.  Then  the  equation  of  motion  becomes 


D^— ^ = P(y)  + F(u))  qZ6  (y-ns) 

djA 


+ G(oo)Z 


3y 


6(y-ns)fe 


We  now  will  introduce  the  Fourier  expansion  for  a 
regular  array  of  6-functions: 


i 


Z6(y-ns)  = (u/2  7r)Ze:*'nuy 

where  u = 2tt/s,  and  we  assume  that  the  net  pressure  and  the 
displacement  have  expansions  of  the  forms 

P - zb  eiPuy+ivy 
p . 

n = za  eiPuy+ivy 

p 

(All  sums  on  p will  run  from  p = -»  to  p = +°°)  . That  is,  we 

assume  that  each  of  P(y)  and  q(y)  is  the  product  of  the  factor 

elvy  and  a function  that  is  periodic  in  y,  with  period  s,  and 

we  expand  this  function  in  a complex  Fourier  series.  The 

coefficients  a and  b depend  on  u and  v.  It  will  be  seen  later 
P P 

that 

, 2 * 2 ,*2. -1/2 
bp  = p ((pu+v)  -k  ) 

where  p*  is  the  water  density,  and  k 
waves  in  the  water. 
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We  define  the  function 

f(w)  = D(w4-k4)-  ca2p*(w2-k*2)"1/2 

and  set: 

Fo  = uF/2  7 r = F/s 
GQ  = uG/2  7r  = G/s 

Then  the  equation  of  motion  becomes: 

Z apf  (pu+v)  e (Pu+V)  iy 

= eivy  + F Z a Z e(nu+Pu+v> ** 

° p P n 

- Gc  Z ap  (pu+v)  Z (pu+nu+v)  e(nu+Pu+v>  ^ 

The  factor  e lv^  can  be  divided  out.  In  the  sums  on  the  right, 
we  make  m = n+p  a new  summation  variable,  in  place  of  n,  and 
change  the  order  of  summation.  Then  we  rename  the  summation 
variable  on  the  left  as  m.  Now  the  coefficients  of  corresponding 
terms  elmu^  can  be  equated: 


a f(mu+v)  = 6 + F Za  - G Za  (mu+v) (pu+v) 

m mo  op  op  ' 'r  ' 


m = -o°,  . . . , +<». 


Now  setting 

A = Z a , P = Z a (pu+v)  , 

p P p P 

we  have 

am  - f’1(v)6mo  + f_1(mu+v)  (FoA-CQ(mu+v)P) 

m = -«>,  . . . , +»  . 
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To  determine  A,  we  sum  this  expression  for  a m over  m;  to  determine 
P,  we  sum  after  multiplying  by  mu+v.  Then  we  obtain 


A = f ” 1 (v)  + UFA  - G U,P 
o o o 1 

P = vf_1(v)  + U^A  - GqU2P 


where  we  have  set 

U.  = Z f ^ (mu+v) (mu+v) * i = 0,  1,  2. 

1 m 

Thus  we  have  obtained  a system  of  two  simultaneous  equations  for 
A and  P: 

(1-FoUq)A  + GqU1P  = f_1(v) 

~ FoUlA  + 0-+GoU2)P  = vf-1(v) 

If  we  define 

T . = Z f \mu+v)m^, 

1 m 

then  the  determinant  A of  this  system  can  be  written  in  the 
alternative  forms 

A . 1 - F0U0  + 0oU2  + FoG0(U^-U0U2) 

' 1 ' FoTo  + 0o(u2T2+2uvT1+v2To)+  F0C0u2(T2-T0T2) 

We  then  have 

A - f"1(v)A‘1(l+G0(u2T2+uvT1)) 

P - f'1(v)A‘1(v^FouT1) 
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The  final  result  is: 

t\  = eivyf_1(v)^l  + ^’1(F0-v2G0+u2F0G0T2)Z  f'1(pu+v)eipuy 


- A'1(uvG0+uF0G0T1)Z  f"1(pu+v)peiPuy | 

It  remains  to  find  P,  and  to  justify  the  form  used  for 
bp.  Consider  a plate  with  water  on  one  side,  but  no  resonators. 
Suppose  that  the  plate  executes  vibrations  with  wave  number  w: 


rj  = eLwy 


at  x = 0.  Then  the  velocity  potential  for  the  acoustic  displace- 
ment at  depth  y is: 

-iw(w2-kV1/2  eiwye-(w2-k*2)1/2x 

and  therefore  the  pressure  at  the  interface  is 

2 *.  2 .*2. -1/2  iwy 
co  p (w  -k  ) e J 

Thus  if  an  arbitrary  displacement  is  Fourier-analyzed  into  terms 
with  wavenumber  v+pu , the  corresponding  pressure  term  will  have 
bp  as  coefficient.  Thus  the  form  taken  for  f(p)  is  correct. 

To  determine  the  displacement,  it  was  not  necessary  to 

specify  whether  the  plate  is  set  into  motion  by  an  incident 

acoustic  wave  in  the  water  or  by  a traveling  wave  of  force, 

distributed  sinusoidally  and  having  wavenumber  v,  applied  to  the 

back  of  the  plate  (the  air  side).  This  distinction  is  essential 

in  determining  the  net  pressure  on  the  plate.  If  a component  of 

displacement  at  wavenumber  w is  due  to  a force  on  the  back,  the 

corresponding  component  of  pressure  is  found  by  multiplying  by 
2 * 2 *2  -1/2 

a)  p (w  -k  ) ' , as  we  have  just  seen.  The  responses  of  the 

regular  array  of  resonators  can  be  analyzed  into  such  traveling 
force  waves,  and  thus  the  effect  on  the  total  pressure  is  known. 
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The  remaining  contribution  to  the  pressure  can  be  found  by  setting 

Fo  = Go  = °*  since  this  contribution  is  independent  of  F and  G . 
Thus  if  we  define  ° ° 


g(w)  = f(w)  J k*2 -w2 , 


we  have 


elVy[l+(to2P*/g<v»+r^{<F0-v2G0+u2FoToT2)S|^) 

o ipuy  I 

+ (uvG  +u  F G T,)Z  £-1  . } 

o o o 1/  g(pu+v)  j 
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APPENDIX  B 

DERIVATION  OF  THE  THEORETICAL  EXPRESSIONS 
V.  Mangulis 


Let  the  incident  plane  wave  pressure  p^  e^0*  in  water 
be  given  by 

Pinc(x’y)  = % exp(-ikx  sin0  - iky  cos0) , (B-I) 

which  is  normalized  in  such  a way  that  the  pressure  magnitude 
would  be  equal  to  one  on  the  plate  itself  if  the  plate  were 
ideally  rigid.  Since  the  structure  repeats  itself  at  inter- 
vals s,  the  total  pressure  p(x,y)  in  water  and  the  velocity 
component  in  the  x-direction  vx(x,y)  should  also  repeat, 
except  for  a phase  factor,  i.e.,  we  should  have  for  integers 
m (suppressing  a time  factor  eiat) 

p(x,y+ms)  = p(x,y)  exp(-ikms  cos0),  (B-2) 

v (x,y+ms)  = v (x,y)  exp(-ikms  cos0) . (B- 3) 

A X 

Consequently  we  can  expand  the  total  pressure  in  a normal 
mode  (or  waveguide  mode)  series: 

P(x»y)  = PinC(x’y) 

'X: 

+ exp(-iky  cos0)  £ An  exp  (inuy  + ,nx)  (B-4) 

n=-« 

where  A^  are  unknown  coefficients  (to  be  determined  by  the 
boundary  conditions  on  the  plate), 


U : 2rr/s, 


^ i (k^-(nu-kcos(J)2]^,  if  0^  ni:-kcos0  ^ k; 


-v,  = 


[ (nu-kcos^)^-k^ , if  | n:-kcos(l  k. 


(B-5) 

(B-6) 
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In  particular,  on  the  plate  itself,  when  x-0, 

00 

p(0,y)  = exp (-iky  cos 0)  £ (Ar+%  6 )exp(inny)  (B-7) 

n=-°° 

where  5 is  the  Kronecker  delta,  6^  = 1 if  n=0,  and  5 « 0 

no  J no 

if  n/0. 

The  velocity  component  in  the  x-direction  vx(x,y) 
and  the  total  pressure  p(x,y)  are  related  by 


vx(x,y)  = (i/po>)  [dp(x,y)/Sx] , (B-8) 

where  p is  the  density  of  water,  and  a>  is  the  radian  frequency. 
At  the  plate 

vx(0,y)  = (i/po>)exp(-ikycos0)  £ ^An'^}no^T|n 


exp  (iniay) 


(B-9) 


which  is  also  the  velocity  of  the  plate  itself  in  the  x-direc- 
tion. 


On  the  plate  the  velocity  v and  the  total  pressure 

x 7 

p must  satisfy  the  differential  equation  (obtained  from  the 
equilibrium  of  forces  and  moments  on  the  plate) 


D 


iu> 


P + 


by 


.hco^v 


(B-10) 


where  is  the  force  per  unit  length,  due  to  a resonator 
(in  the  x-direction),  and  M^,  is  the  moment  per  unit  length 
due  to  a resonator; 

D = Eh3/12 (l-o2)  (B-ll) 


2 b 


„ , 
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where  E is  Young's  modulus  and  a is  Poisson's  ratio;  h is 
the  plate  thickness,  Pp  the  plate  density.  For  the 
resonator  we  let  (note  that  v^/iaj  is  the  plate  displacement 
in  the  x-direction) 


Fl  z (F  vx/ico)  ret  (ms-b,y,ms)  (B-12) 
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I 


rct(s-b,y,s) 


5(y-s+b)-6(y-s), 


(B-16) 


where  5(y)  is  the  Dirac  delta  function. 

If  we  multiply  Eq.(B-15)  by 

exp(-imp.y+ikycos0)  • integrate  over  y from  0 to  s,  and  divide 
by  s,  then  we  obtain  for  m = 0,  +1,  +2,  . ..,  the  infinite 
set  of  equations 


(V*6mo)r>m  [D(n»4-kcos^)4-ppha)2]  - Pa>2(Am+*6mo) 

00 

= £ (An-%Dno)1n(Bnm[(F/b)-(G/b)(n>i-kcosJI)2  (B.17) 

n=-oo 

+ iEnm(G/b)(nu-kcos0) ) 


where 


b/s,  if  n=m, 


Bnm  = < 


(B-18) 


l l-exp[-i(n-m)fiblt  if  n ^ 
2iTi(n-m) 


Enm  = f exP l-i(n“m)ub]-l)/s. 


(B-19) 


We  assume  that  b is  very  small  compared  to  a wave 


length  in  water  or  plate,  and  that  b/s  is  very  small,  so  that 
we  can  simplify  Eqs.(B-17),  (B-18),  and  (B-19)  by  letting 
b -*  0,  which  yields 
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— -i(n-m)n(b/s); 


(B-21) 


and  Eq.  (B-17)  becomes 


(Am-%  &mo)  { r]m[D(mn-kcos)3f)4-pphto2  ] -pa>2 } 

00 

- *“26om  + E <V*6no>V<F/s> 


(B-22) 


- (G/s)  (mp-kcos0)  (ntx-kcos^f)  ] 

While  Eq.(B-17)  consisted  of  an  infinite  set  of 
equations,  and  the  solution  of  that  infinite  set  required  the 
inversion  of  an  infinite  matrix,  by  letting  b-0  we  have  simp- 
lified the  solution  to  the  point  where  we  have  to  solve  only 
a set  of  two  equations,  as  shown  below. 

Note  that  the  following  summations  are  just  some 
constants,  independent  of  m in  Eq.  (B-22): 


(F/s)  £ (V^noK  = bl 


(B-23) 


(G/s)  (An-^nno)Tin(nu-kcos(?)  = b2  (B-24) 


In  terms  of  those  constants  we  can  immediately  write 
down  the  solutions  for  the  unknown  coefficients  Am, 


pu>  +b.  -(mu-kcos0)bo 
A _ko  * om  1 t 

m mo  iim(D(mu-kcosJ?)^-ppha^  ] -pco^ 


(B-25) 


It  remains  to  evaluate  the  two  constants  b^  and 
b2-  We  can  write  down  two  equations  involving  those  two 
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constants  by:  1)  multiplying  both  sides  of  Eq.  (B-25)  by 
(F/s)1^  and  summing  over  m from  -<®  to  which  will  give  b^ 
on  the  left-hand  side;  2)  multiplying  both  sides  by  (G/s)^ 
(mp-kcos0)  and  summing  over  m from  -<»  to  which  will  give 
b2  on  the  left-hand  side.  The  two  equations  can  then  be 
solved  for  b^  and  b2»  which  in  turn  yields  Am  in  Eq.  (B-25), 
which  then  determines  the  pressure  and  velocity  at  the  plate. 

The  final  result  is 


where 

v = kcosf?  (B-28) 

f(nu-v)  = DCnii-v^-Ppha^-pcoV'ln  (B-29) 

g(nu-v)  z ir^fCnu-v)  (B-30) 

Cx  = F/s  - (G/s)v2  + (G/s)u2T2  (B-31) 

C2  = (G/s)u [v+(F/s)uT1 ] (B-32) 


¥ 
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A = l-(F/s)T0  + (G/smAq+ZhvT^h2^) 

-(F/s)(G/s)m.2(T0T2-T12)  (B-33) 

where  for  j * 1,  2,  3 

00 

(-1)j’rj  = Y,  mj  (B-34) 

ms-oo  f(mia-v) 

The  result  derived  here  is  the  complex  conjugate 
of  the  one  derived  in  Appendix  A because  the  time  dependence 
used  here  is  e^CDt>  and  the  one  used  in  Appendix  A is  e"^a)t. 

Equations  (B-26)  and  B(27)  have  been  programmed 
for  a computer.  The  infinite  sums  were  replaced  by  finite 
sums  in  Eqs.  (B-26),  B(27),  and  (B-34);  the  summation  was 
terminated  when  either  the  next  term  in  the  sum  was  less  than 
10'4  times  the  sum  of  previous  terms , or  when  a prescribed 
number  (usually  about  40)  of  terms  had  been  summed.  The 
ratio  of  the  last  term  in  the  sum  to  the  total  sum  is  always 
printed  out  so  that  one  can  always  ascertain  whether  the  sum 
has  converged.  The  convergence  becomes  poorer  as  the  spacing 
s between  resonators  is  increased. 

The  functions  F and  G for  a simple  circular  rod 
resonator  are  ^ 


F = 
G = 


(l/2a)  ESa  tanaL 

EIp  cosgL  sinhgL-sinftL  coshftL 
l+cos0L  coshpL 


(B-35) 

(B-36) 


where  a is  the  radius  and  L the  length  of  the  resonator, 


P4  = co2f)pS/EI 


(B-37) 
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a = ai/(E/p  p)' 


S = Tra 


I = Tra  y4 
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(B-38) 

(B-39) 

(B-40) 


where  we  have  assumed  that  the  resonator  is  made  from  the 
same  material  as  the  plate. 

The  values  which  we  used  in  our  computer  program 
are  (in  MKS  units) 

p - 103  kg/m3,  Pp«  7.7  x 103  kg/m3, 

E z 1.36  x 10"  newtons/m3,  a = 0.29, 

C = 1500  m/sec,  a=  0.06  m, 


L = 0.42  m. 

Other  parameters,  such  as  s,  h,  0r,  and  frequency  f, 
were  allowed  to  vary,  and  are  stated  separately  for  each 
figure  in  the  main  section.  The  values  of  F and  G actually 
used  were  1/10  of  those  given  by  Eqs . (B-35)  and  (B-36)  with 
the  above  values  of  the  parameters  a,  L,  etc. 
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